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212 CYCLIC CODES 



8.2 Matrix Description of Cyclic Codes 

The most elementary description by matrices has been illustrated in the 
previous section. Ug(X) = a r X' + a r _, + • • - + a 0 is the generator 
of the code, then {X-'-g(X)), {X a ~'~ 2 g(X)}, .... (g( X )} are all code 
vectors. Thus, all the rows of the following matrix are code vectors: 



G = 



0 a. 



0 
0 



Or 



0 



0 
a. 



0 



0 
0 

0 

<*oJ 



(8.3) 



Clearly, they are linearly independent, and the rank of G is n - r, which 
is also the dimension of the code. Therefore, by Theorem 2.9 the row 
space of G is the code space. 

As earlier, the following convention will be used. The first k symbols 
coefficients of X"~\ X"~\ .... X"~ k , will be taken as information 
symbols, and the last n-k symbols, coefficients of A"* - * -1 , x n ~ k ~ 2 
... 1, will be taken as parity-check symbols. 

The generator matrix for any cyclic code can be put in modified 
reduced echelon form as follows. Let r ( (X) be the remainder after 
dividing X' by g(X) : 



X i =g(X)g i (X) + r i (X). 



Then 



X s -r i (X)=g(X)q i (X) 

is a code vector. If these polynomials, for i ; = n - I,n-2, ...,n — k 
are taken as rows of the generator matrix, then 

G = [I*,-R], 

where l k is a k x k identity matrix and — R is a k x (n — k) matrix 
whose yth row is the vector of coefficients of -r„_j(X). Then by 
Theorem 3.3 the code is also the null space of the matrix 

M=[R T ,I n _ k ). 

The yth row of H r is the vector of coefficients, of r„_ J .(A*), even for 
jlgn-k. 
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Example. For 
X 3 + X 2 + 1, 



the binary cyclic code generated by f(X) = 



X 6 = g(X)(X 3 + X 2 + X) + X 2 + X 
X s = g(X)(X 2 + X + 1) +X+ I 



X*=g(X)(X+ 1) 

* 3 =<7(*)(l) 
* 2 = 0 (*)(0) 

= <?(*)(0) 

*° = g(X)(0) 



+ X 2 +X+ I 
+ X 2 + 1 

+ x 2 

+ X 

+ 1 



HT 



"1 10 
Oil 

1 1 1 

101 
100 
010 
001 



(8.4) 



and 



G,= 



10001 10 
010001 1 
00101 1 1 
0001 101 



(8.5) 



The matrix G given here is row equivalent to the matrix given in 
Equation 8.1. The code is the same code, not just an equivalent 
code. 

Now consider the dual code that, by Theorem 6.12, is generated 



by(* 7 


- D/9(X) 


= (*- 


IXX 3 + X + l) = X 4 + X 3 + X 2 + l: 


r 6 (X) = 


X 3 + X 2 


+ X 




"1110" 




r 5 (X) = 


X 2 


+ X+ 1 




01 1 1 




U(X) = 


X 3 + X 2 


+ 1 




1101 




r 3 (X) = 


X 3 




H T 2 = 


1000 


(8.6) 


r 2 (X) = 


X 2 






0100 




r x {X) = 




X 




0010 




r 0 (X) = 




1 




_0001 




and 













G 2 = 



1001 1 10 
01001 1 1 
001 1 101 



(8.7) 



The matrices H, and G 2 differ in having the order of columns 
reversed and the rows reversed, and the same is true of G 1 and H 2 . 
The rearrangement of rows does not affect the row space or null 
space. The rearrangement of columns is a result of the fact that two 
polynomials multiply to zero only if the dot product of the vectors, 
with the order of components reversed in one of them, is zero. 
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and by Theorem 6.14, its conjugates are cyclic shifts of this polynomial; 
that is, 

R 2 = a 27 + a 26 + a 25 + a 24 + 0 + a 2 ' + 0 + 0, 

R 27 = (J 9 = 0 + 0 +a 25 + a 24 + a 2i + a 22 + a 2 ' + a 2 ° 

With some effort it can be shown that fi 3 = a 2 ' + cc 2 °. Then the 
conjugates of/? 3 consist of cyclic shifts of the 8-tuple 00 1 0000 1. The 
root f}° can be shown to equal a 27 + a 26 + a 25 + a 2 " + a 21 + a 22 + a 2 * 
+ a 2 ° (See Problem 8.23.) Thus, the matrix H can be rearranged to give 

H = [[/? 9 /? 1 V 15 ^ l6 /3 8 )SW][/? ,0 ^ 5 ^ ,, ^ , W 2 ^ 6 /? 3 ]^ 0 ] 





"01011110" 


'10000100" 


r 




00101 1 1 1 


01000010 


1 




100101 1 1 


00100001 


1 




1 100101 1 


10010000 


1 




1 1 100101 


01001000 


1 




1 1 1 10010 


00100100 


1 




01 1 1 1001 


00010010 


1 




.101 1 1100 


_00001001 


L 



Therefore, the (17,9) code with one information symbol deleted is 
equivalent to a (16, 8) quasi-cyclic code. Because C 2 is nonsingular, the 
matrix [C, C 2 ] can be put in echelon canonical form by multiplying it 
on the left by the inverse of C,. 

It is known that every finite field has a normal basis (Albert 1956, 
p. 1 19). Thus, the generator matrix of every cyclic code can be decom- 
posed into circulants, although many will be trivial. Then, removing 
columns corresponding to elements of order less than n will give the 
generator matrix of a quasi-cyclic code. It may not be possible to put 
such a code in systematic form, however. 

8.15 Codes Based on the Chinese Remainder Theorem 

Let i(X) denote a polynomial of degree k - 1 or less with symbols from 
GF(q m ). Let m^I) denote a polynomial of degree d t with symbols from 
the same field. An interesting class of codes can be constructed from 
the following version of the Chinese Remainder Theorem (Stone 1963). 

Theorem 8.20. The polynomial i(X) can be reconstructed from the 
remainders 



rlX) = i(X) mod m&X) i = 1 , 2, . . . , n, 



(8.56) 
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provided the m { (X) are relatively prime in pairs and that 



(8.57) 



Proof. Any two solutions i x (X) and i 2 (X) to the set of congruences 
(8-56) are congruent modulo m(X) = JJ% % m t (X); for /, (X) - / 2 (;r) is 
divisible by m(X) if and only if it is divisible by each m£X). Since there 
is exactly one polynomial of degree less than k in each residue class 
modulo m{X), there is exactly one solution for the set of congruences 
(8.56), and this must be the correct solution. Q E D 

This theorem can be used to construct a class of random-error- 
correcting codes with symbols from GF(q m ) as follows: Choose as the 
m t (X) the factors of X qm ~ l - I in GF(g m ). Since this polynomial splits 
completely in GF(q m ) 9 

m i (X)=(X-a i ) 9 , = 0,1,. ..,^-2, 

where a is a primitive element of GF(q m ). 

Now calculate the n = q m - 1 residues of the information polynomial 
i(X) which has degree less than k. The set of these polynomials forms an 
(n, k) linear code with symbols from GF(q m ). 

Decoding can be performed as follows: Suppose / errors occur so that 
n -t of the residues are received correctly. Now, i(X) is determined by 

any set of k correct residues. Thus, of the Q ways of determining 

i(X), exactly ^ ^ ^ will agree and give the correct polynomial i(X). It 

is possible that / - 1 transmitted residues have been changed so as to 
be the residues of some other polynomial j(X) of degree k - 1 or less 
which is determined by k - I correct residues and one incorrect 
residue. There are 



determinations which will give the same incorrect answer. 

Now if n - t > k - 1 + /, then the number of correct determinations 
exceeds the number of incorrect ones, so i(X) can be reconstructed 
correctly. That is, if / is chosen such that 




2t + 1 = n - k + 1 
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for a particular code, then ail t errors can be corrected. Since the 
minimum distance is not less than It + 1, and since d ^ n — k + 1 for 
any linear code, it follows that for these codes, 

d = n-k+ 1. 

It turns out that these codes can be put in a cyclic form. The residue 
of i(X) = i k - l X k - i 4- i k . 2 X"- 2 + ..- + /,* + / 0 mod(X - a'), which 
is an element of GF(q m ), is given by 

i(X) = r i +g i (X)(X -a 1 ). 
Thus 

r,- = /W U.. = '•(«•')• 
Let the generator matrix of an (/?, &) code with symbols in GF(q m ) be 
r (a *-i)«-i ( a *-i)»-2 ... a *-« a on 



G = 



(a*- 2 )""' (a*- 2 )"" 2 ••• a*" 2 a 0 



a a or a 



(8.58) 



Multiplying G by the information vector = (i k _ x i k _ 2 " i 0 ) gives the 
/i-tuple 

[ifflk-,-., fWkiW|i-«] 

= [^.,,^.2, ...,r„r 0 ], (8.59) 

which is a code word in the code based on the Chinese Remainder 
Theorem. 

But G is the generator matrix of the cyclic (n, k) code for which 
h(X) has roots oc\ a 2 , . . a*"" 1 . In fact, this is a Reed-Solomon code 
and is discussed at length in the next chapter. 

The decoding procedure presented here for these codes is clearly 
impractical for all but the shortest codes. It requires that the decoder 

perform \ ! J operations for each decoding. Fortunately a more easily 



mechanized method of decoding these very powerful codes has been 
found; it too is discussed in Chapter 9. 

Since a symbol in GF{q m ) can be represented as an /w-tuple over 
GF(q) y these codes are ideally suited to multiple burst-error correction. 
This is discussed in Chapter 1 1. 

It is not necessary to consider linear factors in an extension field; 
polynomials of higher degree over GF(q) also yield interesting codes. 
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Solomon codes and interleaved Reed-Solomon codes, have good burst- 
correcting capabilities. In particular, a Reed-Solomon single-error cor- 
recting code over GF(q b ) interleaved to degree i has exactly the same 
parameters and burst-correcting capability as the block-interleaved 
Burton code (Equation 11.7) with n = q m — I, the maximum possible 
value. (See Problem 11.10.) Note that as a code over GF(q m ) the /-error 
correcting Reed-Solomon code requires 2t check symbols and corrects, 
among other patterns, every burst of length / or less. Thus every Reed- 
Solomon code, as a code over GF(q m ) y is optimal with respect to the 
Reiger bound. It follows, as for the Burton codes, that as codes over 
GF(q m ) they are asymptotically optimal. The multiple-burst-error 
correcting capability of such codes is discussed in Section 11.4. 

Example. A Reed-Solomon code with symbols from GF(2 1 ) and 
/ = 4 could be used to correct all bursts of length 22 or less by 
coding symbols as blocks of 7 binary digits. It would require 
2x4 = 8 blocks or 8 x 7 = 56 binary symbols as parity checks. 
The code length would be 7 x (2 7 - 1) = 889. 

A Fire code to correct all bursts of length 22 or less would 
require 3b — 1 = 65 binary digits as parity checks. The code length 
could be as great as (2 22 — 1) x 43, or approximately 160 million 
binary digits. On the other hand, a burst of length 8 binary digits 
or less could affect at most two elements of GF(2 7 ), and thus the 
Reed-Solomon code could correct any two bursts of length 8 or 
less. It could correct as many as four bursts if each happened to lie 
within a single GF(2 7 ) element. 

From Equation 1 1.7, there exists a (3556, 3500) Burton code with 
b = 22 and z = 12. Note that this code has the same burst-correcting 
ability and number of parity checks as the Reed-Solomon code but 
is four times longer. It cannot correct multiple bursts, however. 

BCH codes can be implemented by the methods described in Chapter 
9. Alternatively, if correction of single bursts only is sufficient, the 
method of Section 11.3 can be used, and requires much less equipment. 
The additional error-correcting ability of the code can be turned to 
error detection (any error pattern that the code would correct, if fully 
utilized, it could certainly still detect) or it can be used to correct bursts 
beyond its guaranteed burst-correcting ability. (See Section 11.3.) 

Theorem 1 1.4. If a code with minimum distance d is used to correct 
multiple bursts of errors, then the set of all burst patterns satisfying 

3d 

l + P^j - 1 (11.8) 
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